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Abstract 

We study the microscopic structure of the Josephson current in a single-mode 
tunnel junction with a wide quasiclassical tunnel barrier. In such a junction 
each Andreev bound state carries a current of magnitude proportional to the 
amplitude of the normal electron transmission through the junction. Tremen- 
dous enhancement of the bound state current is caused by the resonance 
coupling of superconducting bound states at both superconductor-insulator 
interfaces of the junction. The possibility of experimental observation of the 
single bound state current is discussed. 
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The Josephson effect in a tunnel junction deals with coherent transmission of Cooper 
pairs through a tunnel barrier which separates superconducting electrodes. What is the 
mechanism of such a transmission? Conventional theory of the Josephson effect Jl|,0], based 
on the phenomenology of the transfer Hamiltonian model [[|, does not provide adequate 
physical description of this process, e.g. similar to the quantum mechanical picture of 
tunneling of normal electrons. Instead, it treats tunneling as a perturbative transition, 
introducing a matrix element of coupling of electrons in different electrodes proportional to 
the amplitude of single electron tunneling [|J. 

A more realistic description of Josephson tunneling, based on the Bogoliubov-de Gennes 
equation ||, was suggested by Furusaki and Tsukada ||. The crucial role in this picture is 
played by superconducting bound states, similar to Andreev bound states in SNS junctions 
0. A bulk supercurrent, when approaching the tunnel interface, experiences transformation 
into current flowing through superconducting bound states which provide transmission of 
Cooper pairs through the barrier. The bound states are induced in the vicinity of the junction 
by the discontinuity of the superconducting phase, and they appear as a consequence of the 
current In a quantum junction the bound state spectrum consists of a single pair 
of levels per transverse mode, with symmetric position of the levels with respect to the 
chemical potential. The Josephson current is distributed among bound states in such a way 
that each bound state carries a current proportional to the normal electron transparency of 
the junction. 

In this Letter we show that the above picture of quantization of the Josephson coupling 
is valid only for an extremely narrow barrier, and that the picture is qualitatively different 
for any realistic tunnel barrier with large width on an atomic scale. In the latter case, 
the structure of the bound state spectrum is determined by the coupling of superconduct- 
ing surface states situated at the two SI interfaces of the SIS junction. In a symmetric 
junction, the resonance coupling of these surface states provides tremendous enhancement 
of the current flowing through a single bound state, the magnitude being proportional to 
the amplitude rather than the probability of normal electron tunneling. The currents are 
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distributed among the bound states in such a way that they almost cancel each other in 
equilibrium, giving rise to a comparatively small residual current, including the contribution 
from the continuum. This current coincides with the conventional Josephson current given 
by Ambegaokar-Baratoff theory 0. The large current of the single bound state can be re- 
vealed under nonequilibrium conditions when the bound level population is imbalanced by 
means of microwave pumping or tunnel injection. 

We consider for the sake of clarity a single mode quantum constriction with a rectangular 
potential barrier of length L and height V (Fig. 1). The structure is described by the ID 
Bogoliubov-de Gennes equation ||: 
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with the order parameter matrix given by 
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jl = fi — E± is the chemical potential shifted by the energy of the transverse mode, h = 1. 

Consider first an isolated SI surface, assuming in Eq. (1) L = oo. Constructing the 
ansatz for \E\ < A as a superposition of eigenfunctions, with complex wave vectors k± = 



2m(Ji ± i£) (C — V A 2 — E 2 ) in the superconducting region and with decay rates k± 



2m(V — (jl ± E)) in the insulating region, we find the dispersion equation, keeping the 
difference between electron and hole wave vectors to first order, 

E k 2 -\- k 2 ' U 
The surface level in Eq.(3) crucially depends on the finite difference 5k = k + — K- of the 
decay rates of electron and hole wave functions inside the insulator, while dephasing inside 
the superconductor is not important and gives small corrections omitted in Eq. (3). The 
level lies close to the gap edge within the superconducting gap; assuming a typical relation 
among the energies of the problem: A <ti V — jl < jl, one has A — E w (A 3 / '2V 2 )(k / k) 2 . The 
wave function of the surface state decays into superconductor on the characteristic length 



scale I ~ £oA/a/A 2 — E 2 (£ is the coherence length) and into insulator on the characteristic 
length scale I ~ 1/k. 

Proceeding to calculation of the coupling of surface states at finite L < oo in Eq. (1), 
we construct wave functions by means of a transfer matrix formalism f9|, which yields the 
following dispersion relation: 

£ 2 = A 2 cos 2 (^^), 08|±a>0), (4) 

cos a = R + D cos (ft, sin/3 = Dlm(a + a*_), (5) 

where D = ^D + D_, R = ^R + R_, and the transmission D± = |a±|~ 2 and reflection 
R± = 1 — D± coefficients of normal electrons with energy ±E are given in terms of the 
inverse transmission amplitudes 

a± = coshft-i-L — - ( — — — ] sinliK±L. (6) 

2 V k k±J 

Equation (4) is valid for all values of the barrier length L and barrier heights V ||10|| . A 
numerical solution of this equation is presented in Fig. 2, which also shows the corresponding 
wave functions. To get an explicit expression for the level energy we consider barriers which 
are long on the scale of quasiparticle decay but short on the scale of quasiparticle 
dephasing |<5k| _1 : k^ 1 < L <C 1^1 _1 (the last condition is equivalent to (L/£o)(&//c) <C 1). 
In this case, the junction transparency D, the dephasing angle (3 and the coupling factor a 
are all small: D < 1, D + - D_ « -2D6kL, 



a » 2V^V sin W 2 ) + (SkL/2) 2 < 1. (8) 

The last term in Eq. (8) is essential only in the close vicinity of <fi = 0: \<ft\ ~ Since 
the energy dispersion of f3 can be neglected, one finds from Eq. (4): 



E\ = A 2 



sin — 
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(9) 
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where |/?|/2 ± y/D\ sin(0/2)| > 0, D » D, \<f>\ > \Sk\L. 

The solution presented in Eqs. (4), (9) generally consists of two pairs of bands E(<p): 
±E + (4>) and ±E- ((/>), lying inside the superconducting gap symmetrically with respect to 
the chemical potential. The dephasing angle j3 determines the position of the preexisting 
surface states at the single SI interface (cf. Eq. (3)), and the amplitude sfD determines the 
level splitting due to coupling via tunneling. All branches of the bound state spectrum are 



fully developed if \/3\/2 > y/D or kL > In \k/8k\ fl2] . 

The dispersion of the bands E(<f>) determines the current flowing through the bound 
states: / = 2e{dE/d<f>) |l3|,14l . Taking the derivative of Eq. (4), one gets under conditions 
(7), (8) the current of the single level in the form: 

I(E±) = —signE — - — sin0 

x|l ± W 1 V (10) 

2v^D J S m 2 {<p/2) + {5kL/2Y) 



The first term in Eq. (10) is consistent with calculations within the tunnel model, coinciding 
with the Ambegaokar-Baratoff current of a single-mode junction at T = [[J. The second 
term dominates at small transparency VD -C (3/2. It possesses an anomalous square-root 
dependence on the junction transparency and anomalous (close to cos(0/2)) dependence on 
the phase difference, (Fig. 3). This current corresponds to the transition of Cooper pair 
between electrodes with the probability proportional to the amplitude of the normal electron 
transition. From a physical point of view, this enhancement of the tunnel current results 
from resonance coupling of the surface states situated at SI interfaces of the junction, Eq. 
(3). This resonance coupling causes large dispersion of the bound state bands, in analogy 
with the energy level splitting in the Schrodinger symmetric double well potential [IS]. We 
stress that the resonant enhancement concerns only the tunneling of Cooper pairs, while the 
normal electron tunneling is non-resonant. 

Thus, the bound levels in low-transmission tunnel junctions carry a current which tremen- 
dously exceeds the known critical Josephson current. However, this current is not manifested 
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in equilibrium. Indeed, assuming Fermi distribution np{E±) for the occupation numbers of 
the bound levels, one finds with the assumed accuracy the current of all bound states: 

Ibound X; I(E ± )n F (E±) 

signE,± 

= eADtanh(A/2T) sin0. (11) 
This current is twice bigger than the magnitude of the Josephson current calculated for a 



single-mode junction within the tunnel model Q and beyond it [W\. However, in order to 



get the total current one has to take into account the contribution from the continuum: 

I cont = - [ dE[(\t N (<f>)\ 2 - \t A (<P)\ 2 ) - (0 -> -<j>)]n F {E). (12) 

7T J 

The normal transmission amplitude t in Eq. (12) has the form: 

fN _ ~ E[(E - Qa+e^ - (E + flq_ e -*/ 2 ] 

A 2 (cos a + cos j3) - IE 2 cos (3 + 2i£E sin j3 ' 1 ' 

where £ = \JE 2 — A 2 . The Andreev transmission amplitude t A is symmetric in and drops 
out of Eq. (12), which together with Eqs. (7), (8), (13), yields: 

eA 1 
Icont = — —D tanh(A/2T) sin <p = --I bound . (14) 

The calculation presented above, as well as the results of direct numerical calculations in 
Fig. 3, uncover a remarkable fact, namely that the Josephson current in a symmetric junction 
with an extended tunnel barrier results from cancellation of large resonant currents flowing 
through individual bound states. This situation resembles the situation in a normal junction: 
while individual scattering states of electrons carry a finite current, the total current through 
the junction is equal to zero in equilibrium due to cancellation of currents of the modes 
incident from the right and from the left. To reveal the net current of a single mode, one has 
to create a current imbalance, connecting the junction to reservoirs with different chemical 
potentials. In a similar way, creation of current imbalance in the Josephson junction by 
means of nonequilibrium population of the bound states is able to reveal the current of a 
single bound state. For example, it is possible to equalize the level populations within one 
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of the bound level pairs (e.g. ±EL) by means of resonant electromagnetic pumping, as 



suggested in Ref. ||11|| . In principle, one might then suppress the current of this level pair 
and thereby reveal the current of the second pair of levels (±E + ), which will show up in a 
Josephson current enhanced by the factor of \/3\/2\/D, flowing in the same direction as the 
equilibrium current. Suppression of the current of the other level pair ±-E+, by a proper 
choice of the frequency of the pumping field, will show up in an enhanced current flowing in 
the opposite direction. Another possibility is to inject excess quasiparticles into one of the 
bound levels by means of tunnel coupling to an additional normal electrode, i.e. to use a 



three-terminal device similar to the one suggested by van Wees et al. |T7[ in order to reveal 



the currents of individual Andreev states in a SNS junction (see also fll8[ ). 

A crucial condition for observing the current of a single bound state in experiment is 
sufficient resolution of the bound level structure. Since the major part of the bound state 
wave function is concentrated in the electrodes, the main mechanism of level broadening is 
quasiparticle recombination due to electron-phonon scattering in the superconducting banks. 
According to Kaplan et al. []19| this gives an estimate for the level width: T/A ~ at the 
critical temperature, decreasing exponentially at low temperature. The interlevel distance in 
the optimal case of a single-mode junction, according to Eq. (10), is 5E ta A(3\/D, decreasing 
in multimode junctions inversely proportionally to the number of transverse modes. Combi- 
nation of the resolution condition 5E > V with the condition of a long junction, (3 > 2\^D, 
yields for single-mode junction (3 2 > 2T/A. The factor (5 2 depends on the height of the 
tunnel barrier and varies between (A/fj,) 2 < (3 2 < A/ /i in junctions with high (V ~ /i) and 
low (V — fi ~ A) tunnel barriers. In line with these estimates, one might expect to observe 
the effect - giant Josephson current, Eq. (10) - in junctions with a very low tunnel barrier, 
like gate controlled S-2DEG-S devices, or with large ratio A//i, like high T c junctions. 

A more realistic approach, however, would be to increase the energy level splittings, 
increasing the j3 factor by making a SNI well at the interface |[20|| . For a symmetric SNINS 



structure we find the same result for the bound state spectrum as Eq.(4), but with a factor 
(3 = 4d/£o, d -C £cb where d is the width of normal region, corresponding to preexisting 



surface state in SNI well pTJ. This can easily provide a sufficiently large /5-factor, say ft ~ 0.1 



0- 

In conclusion, we have studied a microscopic mechanism of Josephson current transport 
in single-mode tunnel junctions with wide quasiclassical tunnel barriers. We found that 
each Andreev bound state in such junctions carries a giant current proportional to the 
amplitude of normal electron transmission through the junction. Experimental observation 
of the current of a single bound state is possible under nonequilibrium conditions. 
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stages of this work. This work has been supported by the Swedish Natural Science Research 
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FIGURES 

FIG. 1. Single-transverse-mode adiabatic junction modeled as a 1-dim SIS junction with a wide 
barrier of length L and height V. 

FIG. 2. Andreev level energies and wave functions in a single-mode tunnel junction. Left part: 
Andreev level energies E{4>) > as functions of phase difference (f> across the junction. Full lines: 
Long SIS junction (L = 1, V = 50); Dashed line: Short SIS junction (L = 0.2, V = 270) having the 
same transparancy D as the long junction but showing only a single Andreev level, E + [12]. Right 
part: Andreev wave functions |^| 2 = |it| 2 + \v\ for the long SIS junction at = 7r; each of the u 
and v components are symmetric or antisymmetric around x = 0, making E± nearly "bonding" 
or " antibonding" states. In this figure the choice of parameters is beyond the approximate Eqs. 
(8), (9) to achieve suitable resolution. 

FIG. 3. DC Josephson currents in a single-mode tunnel junction at T = (long SIS junction: 
D « 2.4 x 10~ 9 , (3 ss 3 x 10~ 2 ). i+ = I(—E+) and I- = I(—E_) are giant Josephson currents 
associated with individual Andreev levels, hound = 1+ + I- is the small residual current from 
compensating giant bound state currents. I CO nt is the continuum contribution to the current (Eq. 
(12)). I j = hound + hont is the (total) Josephson current of the junction. 
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